Terahertz pulse generation via optical rectification in photonic crystal microcavities 
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Using a 3D fuUy-vectorial nonlinear time-domain analysis we numerically investigate the gen- 
eration of terahertz radiation by pumping a photonic crystal microcavity out of resonance. High 
quality factors and a quadratic susceptibility lead to few-cycle terahertz pulses via optical rectifica- 
tion. Material dispersion as well as linear and nonlinear anisotropy is fully accounted for. 



In the last decades the increasing need of terahertz 
(THz) sources for biomedical applications and for spec- 
troscopy has dragged the researchers' attention toward 
new methods to generate radiation in this wavelength 
region. Several ideas have been implemented, from 
semiconductor-based fotoconductive antennas to optical 
rectification in crystals. _2, 3, 4, 5j Resonant structures 
have been proposed in order to improve THz genera- 
tion efficiency; in this context Photonic Crystals (PC) 
offer various degrees of freedom for device design and 
optimization, both in terms of dispersive features and 
nanoscale integration. PC consist of a high index- 

contrast periodic distribution in either one, two or three 
spatial dimensions: they severely affect light propagation 
when its wavelength is comparable to their characteristic 
period and, since electromagnetic radiation can be con- 
fined in extremely small volumes for long times, PC lend 
themselves to a new class of microcavities with engineer- 
able Quality Factors (Q) as high as 10^. [11113 Several 
authors have recently recognized the advantages encom- 
passed by the use of both linear and nonlinear resonant 
structures and PC in the THz region. [11 El El [Hi 
In this Letter we numerically show that light modula- 
tion at terahertz frequencies can be achieved by resorting 
to a high-Q PC microcavity fed by an out-of-resonancc 
pump. Exciting a long-lived PC cavity-mode with a sin- 
gle source slightly detuned from resonance, taking ad- 
vantage of transient oscillations in a quadratic medium 
we obtain THz radiation via optical rectification. The 
dynamics behavior depends on both detuning and Q. 
In order to simulate the generation process we developed 
a three dimensional parallel Finite Difference Time Do- 
main (FDTD) code, able to account for second-order non- 
linearities as well as material anisotropy and dispersion. 
To this extent Maxwell equations are coupled to an os- 
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cillating dipole with an anharmonic term: |16| 



-D : PP = eo(es - l)L^jiE. 



(1) 

The last equation of |^ assigns the material proper- 
ties: P is the polarization due to a single-pole Lorentz- 
dispersion centered in luq, and the other parameters yield 
a refractive index n ~ 3.5 at A = 1.55^m, with static 
permittivity = 11.7045, loss coefficient 70 = 3.8e8, 
WR = 1.1406el6 and = 1.0995el6 in MKS units. The 
tensor D provides the quadratic nonlinear response, and 
includes a 7r/4 rotation of the crystallographic axes (with 
respect to the spatial coordinates) in order to generate a 
second-order polarization from a linearly polarized input. 
All linear and nonlinear prop erties, including dispersion, 
match those of AlGaAs. (Iq . The microcavity consists 
of a void in a six-hole PC-wire (w = 450 nm, h — 270 
nm and I = 4/Ltm, see Fig. 1 inset) designed to yield 
single-mode TE-like propagation in both input and out- 
put waveguides. [l3| For a — 450 nm and d = 0.35a, 
a resonant state appears near the middle band-gap at 
Ao = 1.447/im, with Q=373. Fig.l graphs the oscillation 
spectrum, obtained by exciting the device with a single- 
cycle (broad-band) pulse and Fourier-transforming the 
field as probed inside the structure. The arrow indicates 
the off-resonance input wavelength. 
To reveal the transient due to microcavity excitation, we 
launched (in the input waveguide) an "mnm" pulse, E3 
which ensures precise control of both leading and trailing 
edges of the pulse and its bandwidth, i.e. avoids spurious 
spectral artifacts. 

To theoretically assess the numerical results, we also em- 
ployed coupled mode theory in the time domain for the 
case of a pumped cavity with feeding waveguide. E3 The 
decay rate l/r = 2/(tJoQ) = l/^o + 2/re includes both 
internal losses (I/tq) and those (I/te) due to the external 
load, i.e. input and output waveguides. For a step-like 
excitation s{t) = for t < and s{t) = sexp[jujt] for 
t > 0, the time evolution of the cavity mode a{t) is given 
by (the input power is Pinput = |sp): 



a{t) 



j{uj - luq) + l/r 



pj^t _ pU^O-i/r}t-\ 



(2) 



2 




1.2 1.4 1.6 

Wavelength (ji m) 



FIG. 1: Spectrum of the electric field (y component) probed 
inside the cavity. The arrow indicates the out-of-resonance 
input wavelength. The inset is a sketch of the device geome- 
try. 



Fig. 12 displays the time evolution of the y-component 
of the electric field inside the cavity for two different 
wavelengths. The top panel corresponds to an excita- 
tion at AA = \input — Ao = 15nm, the bottom one at 
A A = 20nm. As expected from Eq. ^ the larger the 
shift from resonance, the higher the frequency of the os- 
cillating transient. THz waves result from optical rec- 
tification of such transient inside the cavity, due to the 
locally pronounced radiation intensity. Such generation 
can be written in term of the quadratically nonlinear po- 
larization = x'^^EkEi (using contracted notation), 
with x^^^ the pertinent component of the susceptibility 
tensor, mixing intra-cavity fields and Ei, correspond- 
ing to the spatial components /fe(r) and /z(r)of the defect 
mode, respectively. Naming 5 = ojq — u the frequency 
generated via the interaction, the resulting amplitude of 
the THz nonlinear polarization is: 



^THz 



To/Te 



ijiput 



(1 + 2TjTe)2 + ,5^2 



(3) 



The temporal dynamics is governed by the cavity 
decay-time, and provides few-cycle pulses at THz. Since 
the latter wavelengths are in the millimeter range, the ra- 
diation is unconfined and barely "senses" the structure. 
If required, however, it could be trapped e.g. by a hol- 
low pipe waveguide. "20] Fig. O displays the generated 
THz power evaluated by FDTD from the power density 
spectra (PDS) and at two detuning 5. The curves clearly 
exhibit a quadratic dependence from the input power at 

the lines being quadratic fits. 
While the final conversion efficiency depends on the 
(eventually) adopted THz guiding geometry, the case 
A = 1.462/im (corresponding to ^ = 2.\THz) is clearly 
more efficient than for A = 1.467/im [5 = 2.8THz), as 



expected (see Eq. 

To demonstrate how this scheme works for a pulsed ex- 
citation, in our FDTD code we fed the input waveguide 
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FIG. 2: Electric field evolution of the signal inside the cavity 
when pumped by a step-like excitation at A = 1.462^m and 
A — 1.467/im, corresponding respectively to AA = 15 nm and 
AA = 20 nm. 



with two "mnm" pulses (m — 6, n = 600), mimicking 
a Return-to-Zero modulation of the optical pump. The 
top panel of Fig. Qlshows the y component of the electric 
field inside the cavity when pumped at A = lA62fim. 
To extract the THz portion of the electromagnetic field 
we performed the spectrogram of the TM component, al- 
lowing it to follow the time evolution of the signal at the 
desired frequency. The bottom panel of Fig. 21 plots the 
THz pulses obtained through this procedure. 
Finally, to quantify the generation efficiency 77, we eval- 
uated the intensity of both optical and THz signals, in- 
tegrating the PDS of TE and TM components around 
ujQ and S, respectively. This yielded a promising rj — 
2.3 X 10-^ for S = 2.1THZ and 77 = 1.3 x lO"'^ for 
6 — 2.8THz, respectively. Note that, in this unopti- 
mized 3-hole PC wire, Q — 373 is far below the values 
rece ntly obtained in PC cavities, as high as Q = 1 x 10^. 

mm 

The best performance can be achieved when max- 
imizes the last factor in Eq. Q, hence THz generation 
becomes linearly dependent on Qo and the final efficiency 
grows quadratically to values competitive with those re- 
cently achieved. Q 

In conclusion, we numerically investigated THz sources 
based on optical rectification of the transient response of 
a PC-wire microcavity. This optically pumped device, 
amenable to be incorporated in THz circuitry, can be 
tailored to specific requirements in terms of frequency 
oscillation and efficiency. 



3 




FIG. 4: Two mnm pulses at A = 1.462/im (top panel) giving 
rise to two THz pulses after filtering (bottom panel). 
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